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ABSTRACT 



The high sensitivity of the Fermi-LAT (Large Area Telescope) offers the first 
opportunity to study faint and extended GeV sources such as pulsar wind nebulae 
(PWNe). After one year of observation the LAT detected and identified three pul- 
sar wind nebulae: the Crab Nebula, Vela-X and the PWN inside MSH 15-52. In 
the meantime, the list of LAT detected pulsars increased steadily. These pulsars are 
characterized by high energy loss rates (E) from ~3 x 10 33 erg s _1 to 5 x 10 38 erg 
s _1 and are therefore likely to power a PWN. This paper summarizes the search for 
PWNe in the off -pulse windows of 54 LAT-detected pulsars using 16 months of sur- 
vey observations. Ten sources show significant emission, seven of these likely being 
of magneto spheric origin. The detection of significant emission in the off-pulse in- 
terval offers new constraints on the 7-ray emitting regions in pulsar magnetospheres. 
The three other sources with significant emission are the Crab Nebula, Vela-X and a 
new pulsar wind nebula candidate associated with the LAT pulsar PSR J1023— 5746, 
coincident with the TeV source HESS J1023-575. We further explore the association 
between the H.E.S.S. and the Fermi source by modeling its spectral energy distribu- 
tion. Flux upper limits derived for the 44 remaining sources are used to provide new 
constraints on famous PWNe that have been detected at keV and/or TeV energies. 

Subject headings: catalogs - gamma rays: observations - pulsars: general 
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1. Introduction 



Since the launch of the Fermi Gamma-Ray Space Telescope (formerly GLAST) the number of 
detected pulsars in the gamma-ray domain has dramatically increased. The list of LAT pulsars now 
contains 56 bright sources and certainly many more will be detected in the coming months. Yet 
most of the pulsar spin-down luminosity is n ot observed as pulsed phot on emission and is instead 
carried away as a magnetized particle wind ( Gaensler and Slane 2006|) . The deceleration of the 
pulsar-driven wind as it sweeps up ejecta from the supernova explosion generates a termination 
shock at which the particles are pitch-angle scattered and further accelerated to ultra-relativistic 
energies. The pulsar wind nebula (PWN) emission, including synchrotron and inverse Compton 
components, extends across the electromagnetic spectrum from radio to TeV energies. PWNe 
studies can supply information on particle acceleration mechanisms at relativistic shocks, on the 
evolution of the pulsar spin-down and, at later phases, on the ambient interstellar gas. 

Despite the detection of 271 sources, EGRET could not firmly identify any PWNe besides 
the bright Crab Nebula. Most of the 170 unidentified EGRET sources at low Galactic latitudes 
(|6| < 5°) are associated with star-forming regions and hence may be pulsars, PWNe, supernova 
remnants (SNRs), winds from massive stars, or high-m ass X-ray binaries (e.g. Kaaret & Cottam 



1996 



Yadigaroglu & Romanilll997l : lRomero et all 1999). The early LAT observations (|Abdo et al 



2010a|) show that Fermi is detecting many nearby young pulsars. All Fermi-LAI pulsars have a 
high energy loss rate (E), ranging from ~3 x 10 33 erg s _1 to 5 x 10 38 erg s _1 . About a third of 
these pulsars are associated with pulsar wind nebulae candidates observed in the TeV energy range 
by Cherenkov telescopes. These pulsars are thus likely to power a PWN detectable by Fermi. 
However, up to ~10 GeV, the pulsed emis sion dominates the signal from the associated PWN, as 
can be seen with the example of Vela-X (|Abdo et alj|2010cl) . A search for pulsar wind nebulae 
candidates around all detected Fermi-LAI pulsars thus requires that one first removes the pulsar 
signal, thereby selecting only the unpulsed photons. 

Here we report on the analysis of the off -pulse emission of 54 pulsars detected in the gamma- 
ray do main by Fermi-LXY using 16 mon ths of survey observations : 45 pulsars^ reported in lAbdo et al 



(I2010a|) . the 8 new blind search p ulsars (ISaz Parkinson et aU2010 ) and the millisecond pulsar PSR 



J0034- 0534 (lAbdo et alJl2010eh . The study of the PWN in MSH 15-52, reported in lAbdoetal. 



(|2010h|) . did not require the selection of off-pulse photons. Therefore, its associated pulsar PSR B 1509— 58 
is not added to our list of sources. 



The primary objective of this study is to examine the properties of the off -pulse emission of 



'The pulsar PSR J1747— 2958 and its associated off-pulse emission will be studied individually due to its proximity 
to the Galactic center 
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each pulsar and attempt to detect the potential emission associated with its pulsar wind nebula. 
This first population study in high energy gamma-rays allows us to address astrophysical questions 
such as: 

• Do we see pulsar wind nebulae in all Fermi-LAT gamma-ray pulsars ? If not, is it because 
of some specific properties of the pulsar wind or of the ambient medium ? 

• What is the gamma-ray efficiency of PWNe and what physical parameters determine its 
value in addition to the spin-down luminosity of the pulsar ? 

• What fraction of TeV PWNe candidates are detected in the Fermi-LAT energy range ? 

The structure of the paper is as follows: section [2] describes the LAT, sections [3] and |4]present 
the timing and spectral analyses, while the results are described in section |5J Finally, our conclu- 
sions are summarized in section [6J 



2. LAT description and observations 

The LAT is a gamma-ray telescope that detects photons by conversion into electron-positron 
pairs and operates in the energy range between 20 MeV and 300 GeV. It is made of a high- 
resolution converter tracker (direction measurement of the incident gamma-rays), a CsI(Tl) crystal 
calorimeter (ener gy measurement) an d an anti-coincidence detector to identify the background of 
charged particles (|Atwood et al]|2009|) . In comparison to EGRET, the LAT has a larger effective 
area (~ 8000 cm 2 on-axis above 1 GeV), a broader field of view (~ 2.4 sr) and superior angular 
resolution (~ 0.6° 68% containment at 1 GeV for events converting in the front section of the 



tracker). Details of the instrumen ts and data process ing are given in 
on-orbit calibration is described in lAbdo et al . (2009a). 



Atwood etalJ (120091) . The 



The following analysis used 16 months of data collected from August 4, 2008 (MJD 54682), 
to December 16, 2009 (MJD 55181), except for some pulsars for which portions of the observation 
period were rejected due to inadequate pulsar ephemerides, reported in Table [U The Diffuse class 
events were selected (with the tightest background rejection). From this sample, we excluded 
gamma-rays with a zenith angle larger than 105° because of the possible contamination from Earth 
limb photons. We used P6_V3 post-launch instrument response functions (IRFs) that take into 
account pile-up and accidental coincidence effects in the detector subsystems B 



2 See http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/CiceronelATJRFs/IRF_overview.html 
for more details 
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3. Timing analysis 



Most of the pulsars detected by Fermi-LAI are bright point so urces in the gamm a-ray sky up 
to ~10 GeV, though the Vela pulsar is well detected up to 25 GeV (Ubdo et alJl2010bh . The study 
of their associated pulsar wind nebulae thus requires us to assign phases to the gamma-ray photons 
and select only those in an off -pulse window, thereby minimizing contributions from pulsars. We 
phase-folded photon dates using both the Fermi plug-in provided by the LAT team and distributed 
with the TEMP02 pulsar timing packaged, as well as accurate timing solutions either based on 



radio timing o bservations made at th e Jodrell Bank (IHobbs et al] 120041), Nancay (ITheureau et al 



2005), Parkes (IWeltevrede et a 



20091) or Green Bank telescopes (IKaplan et al. or on gamma 



ray data recorded by the LAT (|Ray et al.ll2010h . Whenever possible, data from multiple radio tele- 



scopes were combined to build timing solutions, thereby improving their accuracy and expanding 
their time coverage. 

The origins of the timing solutions used in this analysis can be found in Table CD For each 
pulsar, we list the observatories that provided the data used to build the timing model. For some 
pulsars, we could not produce a timing solution providing accurate knowledge of the rotational 
phase over the whole observation range due to glitch activity. In these cases the time intervals 
over which we lost phase-coherence were rejected. These intervals are given in the last column . 
Also listed in Table [D are the pulsar distance (see lAbdo et al.l l2010al : ISaz Parkinson et al.l I201U : 



Theureau et al.ll2010L for PSR J0248+6021) and the definition of the off -pulse region. These off- 



pulse i nt ervals are chosen using th e definition reported in previous Fermi-LAI studies (lAbdo et al. 



2010a e; Saz Parkinson et aDl201o|) but narrowed slightly to minimize the contamination by pulsed 



photons. A few notes on these timing solutions: 



• The rms of the timing residuals is below 0.5% of the pulsar's rotational period in most cases, 
but ranges as high as 3.6% for PSR J 1846+09 19 which has one of the lowest gamma-ray 
fluxes. This is adequate for the analysis performed for this paper, as timing solutions are 
used only for rejecting pulsed photons. 

• Glitch activity was observed for 12 pulsars over the time range considered here. These 
pulsars are labeled with a g in Table [TJ In all cases it was possible to model the glitch 
parameters in such a way that all the timing data could be used except for PSRs J0205+6449, 
J1413— 6205 and J1813— 1246 where some data had to be rejected as shown in Tabled] 

• Timing solutions were built using radio timing data for all radio-emitting pulsars except 
PSRs Jl 124-5916, J 174 1-2054, J 1907+0602 and J2032+4127. The first is very faint in 



3 http://sourceforge.net/projects/tempo2/ 
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radio a nd was more easily timed in gamma -rays. The three others were discovered re- 
cently (|Camilo et al]|2009l : lAbdo et al.ll2010gh and radio timing observations were therefore 
unavailable for most of the gamma-ray data considered here. For pulsars without radio emis- 
sion, timing solutions were built using the data recorded by the Fermi-hAT only. 



4. Analysis of the Fermi-hAT data 



The spectral analysis was performed using a maximum-likelihood method (M attox et al. 1 1996k 



implemented in the Fermi Science Support Center science tools as the "gtlike" code. This tool fits a 
source model to the data along with models for the diffuse backgrounds. Owing to uncertainties in 
the instrument performance still under investigation at low energies, only events in the 100 MeV - 
100 GeV energy band are analysed. We used the map cube file gllJem_v02.fit to model the Galac- 
tic diffuse emission together with the corresponding tabulated model isotropicJem_v02.txt for the 
extragalactic diffuse and the residual instrument emission^. The off-pulse spectra were fit with a 
power-law model assuming a point- sour ce located a t the position of the pulsar. Nearby sources 



in the field of view are extracted from lAbdo et al.l (I2010i|) and taken into account i n the study. 



Sourc es within 5° of the pulsar of interest and showing a significant curvature index (|Abdo et al 



2010i|) were left free for the analysis assuming an exponential cut-off power-law model, while 
other neighbouring sources were assigned fixed power-law spectra unless the residuals showed 
clear indication of variability from the 1FGL catalog. 

To provide better estimates of the source spectrum and search for the best PWN candidates, 
we split the energy range into three bands, from 100 MeV to 1 GeV, 1 to 10 GeV and 10 to 100 
GeV. The uncertainties on the parameters were estimated using the quadratic development of the 
log(likelihood) around the best fit. In addition to the spectral index T, which is a free parameter in 
the fit, the important physical quantities are the photon flux F .i_ino (in units of ph cm~ 2 s _1 ) and 
the energy flux Go.i-ioo (in units of erg cm~ 2 s _1 ). 

,100 GeV dN 

■Fo.i-ioa = / lTT djE ' and W 

J 0.1 GeV Uh/ 

,100 GeV dN 

Go.i-ioo = / E—dE. & 

J 0.1 GeV Q-Hs 

These derived quantities are obtained from the primary fit parameters and corrected for the de- 
creased exposure represented by the restriction to the off-pulse phase window. Their statistical 
uncertainties are obtained using their derivatives with respect to the primary parameters and the 



Available from http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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covariance matrix obtained from the fitting process. The estimate from the sum of the three bands 
is on average within 30% of the flux obtained for the global power-law fit. 

An additional difficulty with this search is that we must address cases where the source flux 
is not significant in one or all energy bands. For each off-pulse source analysed, gtlike provides 
the Test Statistic TS = 2Alog(likelihood) between models with and without the source. The TS is 
therefore a measure of the source significance, with TS = 25 corresponding to a significance of just 
over 4.5a. Many sources have a TS value smaller than 25 in several bands or even in the complete 
energy interval. In such cases, we replace the flux value from the likelihood analysis by a 95% 
C.L. upper l imit in Tables [2] and [3l These upper limits were obtained using the Bayesian method 



proposed by lHelend (I1983|) . assuming a photon index T = 2. 



All fluxes and upper limits as well as the statistical uncertainties obtained using this procedure 
are summarized in Tables [2] and [3] and were all cross-checked using an analysis tool developed by 
the LAT team called "Sourcelike". In this method, likelihood fitting is iterated through the data 
set to simultaneously optimize the position and potential extension of a source, assuming spatially 
extended source models and taking into account nearby sources as well as Galactic diffuse and 
isotropic components in the fits. The results from this analysis, assuming a point-source model, 
are consistent with those from the likelihood analysis. 

In addition to this cross-check using sourcelike, we performed a second fit to the data with 
gtlike incorporating the results from the first maximum likelihood analysis for all sources other 
than the one being considered, so it has a good representation of the surroundings of the source. 
This step returns a full Test Statistic map around each source of interest. These TS maps do not 
show any extended emission that could contaminate our source of interest (due to badly resolved 
diffuse background) at a TS level higher than 16. 



5. Results 

Pulsar wind nebulae candidates were selected using two different criteria: 

1. TS > 25 in the whole energy range (100 MeV - 100 GeV) 

2. TS > 25 in one of the three energy bands (100 MeV - 1 GeV, 1 - 10 GeV, 10 - 100 GeV) 
As can be seen from Tables [2] and [31 10 of the 54 pulsars studied here sati sfy one of these de 



tection criteria: J0034-0534, J0534+2200 associated with the Crab Nebula (lAbdo et all 201 Od) 



J0633+ 1746 (Geminga), J0835-4510 associated with the Vela-X pulsar wind nebula (lAbdo et al 



2010ch . J1023-5746, J1813-1246, J1836+5925, J2021+4026, J2055+2539 and J2124-3358. A 
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detailed study of the Cr ab Nebula with a m odel adapted to the synchrotron component at low en- 
ergy was performed in lAbdo et al.l (I2010d|) and enabled its clear detection and identification by 
Fermi-LAT. Similarly, a detailed morphological and spe ctral analysis allo wed the detection of the 
extended emission from the Vela-X pulsar wind nebula (lAbdo et al.ll2010q) . 



Aside from the Crab and Vela pulsars, J1023— 5746 is the only candidate that shows off-pulse 
emission predominantly above 10 GeV, whereas the 7 others are mainly detected at low energy 
(below 10 GeV) which suggests a low energy cutoff and therefore a pulsar origin. To provide 
further details on these 7 sources and ensure that the emission detected in the off-pulse interval 
does not have a pulsar origin, we re-fitted all candidates using an exponential cutoff power-law 
spectral model; the results on the off-pulse emission of J1023— 5746 are presented in section I5T21 



5.1. Magnetospheric emission in the off-pulse window 

We explored whether the exponential cutoff power-law spectral model is preferred over a 
simple power-law model by computing TS cuto ff = 2Alog(likelihood) (comparable to a x 2 distri- 
bution with one degree of freedom) between the models with and without the cutoff. The pulsars 
J0633+1746, J1836+5925, J2021+4026 and J2055+2539 present a significant cutoff (TS cuto s > 
9), J2 124— 3358 being at the edge. Pulsars with TS cn to$ < 9 have poorly measured cutoff en- 
ergies; in this case (for J1813 — 1246), we report in Table|4]the fit parameters assuming a simple 
power-law. We also determined if an extended uniform disk model (compared to the point-source 
hypothesis) better fits the data for each candidate. For this step, we used sourcelike and computed 
TS ext = TS 'disk — TS 'point- We did not find any candidates with significant extension (TS ext > 9). 

The Fermi-LAT spectral points for each source listed in Table|4]were obtained by dividing the 
100 MeV - 60 GeV range into 6 logarithmically-spaced energy bins and performing a maximum 
likelihood spectral analysis in each interval, assuming a power-law shape for the source with a 
fixed photon index. The results, renormalized to the total phase interval, are presented in Figures [T] 
and|2]together with the maximum likelihood fit in the whole energy range, assuming an exponential 
cutoff power-law (dashed blue line) or a power-law (dotted-dashed green line). This analysis is 
more reliable than a direct fit to the spectral points of Figured] and [2] since it accounts for Poisson 
statistics of the data. 

Three different systematic uncertainties can affect the results derived with this analysis. The 
main systematic at low energy is due to the uncertainty in the Galactic diffuse emission. Different 
versions of the Galactic diffuse emission generated by GALPROP were used to es timate this error 



in the case of the supernova remnants W51C and W49 (lAbdo et al.ll2009bll2010k|) . The difference 



with the best fit diffuse model is found to be < 6%. Therefore, we estimated this systematic 
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enor by changing the normalization of the Galactic diffuse model artificially by ±6%. The second 
uncertainty, common to every source analyzed with the LAT, is due to the uncertainties in the 
effective area. This systematic is estimated by using modified instrument response functions (IRFs) 
whose effective area bracket that of our nominal IRF. These 'biased' IRFs are defined by envelopes 
above and below the nominal dependence of the effective area with energy by linearly connecting 
differences of (10%, 5%, 20%) at log(E) equal to (2, 2.75, 4), respectively. The third systematic 
is related to the morphology and spectrum of the source. Taking a power-law spectral shape and 
a point-source morphology at the pulsar position are strong assumptions that can affect the flux 
and the spectral indices of the off-pulse compon ent derived with th is simple analysis, as has been 
demonstrated for the case of the Vela-X pulsar (Abdo et al.ll2010q) . A more detailed analysis of 
each source is beyond the scope of this paper and must be handled on a case by case basis. We 
combine the other two systematic errors in quadrature to estimate the total systematic error at each 
energy and propagate it through to the fit model parameters reported in Table HI 

The lack of extended emission and the significant spectral cutoffs at low energies (from 0.43 
to 1.71 GeV) suggest that the off-pulse emission detected by Fermi-LAT is lkely magneto spheric 
and that we do not observe pulsar wind nebulae for J0633+1736, J1836+5925, J2021+4026, 
J2055+2539 and J2124— 3358. This was already su ggested in previous Fermi-LAI publica tions 
on the first two pulsars, J0633+1746 (kbdo et alJboiOjh and J1836+5925 (Abdo et al.lboiOfh . 



The cases of the Fermi-LAT pulsar PSR Jl 8 13 - 1246 and the millisecond pulsar J0034-0534 
are harder to handle due to the limited statistics. For J0034— 0534 an unpulsed component of 
emission from particle acceleration in the wind termination shock might be expected since this 
pulsar is in a binary system, though the two-pole caustic model also predicts a faint signal in 
the off-pulse window of J0034— 0534. In the case of Jl 8 13-1246, which shows a steep spectrum 
with no significant cutoff, we cannot rule out the PWN origin with the current statistics. We 
therefore cannot definitely determine the origin of the emission detected by Fermi-LAT for these 
two candidates. 



5.2. A plausible pulsar wind nebula candidate powered by PSR J1023— 5746 

5.2.1. Fermi-LAr results on the off-pulse emission of PSR J1023—5746 
In 2007, H.E.S.S. reported the detection of very high-energy gamma-rays fr om an extended 



sourc e, HESS J 1023— 575, in the direction of the young stellar cluster Westerlund 2 (|Aharonian et al 



2007|) . Four scenarios to explain the TeV emission were suggested: colliding stellar winds in the 
WR 20a binary system (although this scenario can hardly reproduce the observed source exten- 
sion of 0.18°), collective effects of stellar winds in the Westerlund 2 cluster (although the cluster 



-12- 



angular extent is smaller than that of the very high energy gamma-ray emission), diffusive shock 
acceleration in the wind-blown bubble itself, and supersonic winds breaking out into the interstel- 
lar medium. Recently, Fermi-LAI discovered the very young (characteristic age of 4.6 kyr) and 
energetic (spin-down pow er of 1.1 x 10 37 erg s" 1 ) pulsar J 1023-5746, coincident with the TeV 



source HESS J 1023 -575 ( Saz Parkinson et al.l l2010). 



As noted above, J1023— 5746 is the only candidate that does not show any off-pulse emission 
below 10 GeV, whereas its signal above 10 GeV is > 3a. Therefore, an exponential cutoff power- 
law model, as used for the 7 other candidates, will not represent the data properly. For these 
reasons, we decided to analyze this source separately. 

We searched for a significant source extension using sourcelike with a uniform disk hypothesis 
(compared to the point-source hypothesis). The difference in TS between the uniform disk and 
the point-source hypothesis is negligible which demonstrates that the two models fit equally well 
with the current limited statistics. We have also examined the correspondence of the gamma-ray 
emission with different source shapes by using gtlike with assumed multi-frequency templates. For 
this exercise we compared the TS values of the point source, uniform disk and Gaussian spatial 
models with values derived when u sing a morphological template from the H.E.S.S. gamma-ray 



excess map (lAharonian et al.H2007|) . We did not find any significant improvement (difference in 
TS ~ 3) between the different models and we therefore cannot rule out a simple point source 
morphology. 

To further investigate the off-pulse spectrum and avoid reliance on a given spectral shape, we 
derived the spectral points by dividing the 100 MeV - 100 GeV range into 6 logarithmically-spaced 
energy bins and performing a maximum likelihood spectral analysis in each interval assuming a 
point source at the position of the pulsar (as explained in section IBTTI) . The result, renormalized to 
the total phase interval, is presented in Figure [3] with a red point and arrows. The signal is only 
significant above 10 GeV and is consistent with the H.E.S.S. spectral points. 



5.2.2. Broad-band modeling 

The connection between the GeV flux as observed by Fermi and the TeV flux as seen by 
H.E.S.S. supports a common origin for the gamma-ray emission. The extension of the H.E.S.S. 
source, the off-pulse Fermi signal, and the energetics of this young pulsar point towards a pulsar 
wind nebula origin. The very large number of PWNe detected in the TeV energy range (the most 
numerous class of Galactic TeV sources) and the significant number of PWNe associated with 
Fermi-LAI pulsars make this scenario highly probable. 

Analysis of CO emission and 21 cm absorption along the line of sight to Westerlund 2 gives a 
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kinematic distance of 6.0 ± 1.0 kpc to the star cluster ( DameH2007|) . The assumption that TeV 
emission stems from the pulsar associating PSR J 1023— 5746 with Westerlund 2 is problematic, 
however. The 8' separation of the pulsar and the cluster imply an extremely high transverse velocity 
of ~ 3000 km s" 1 for a 6 kpc distance and the pulsar's characteristic age of 4.6 kyr. In addition, the 
0.18° extension of HESS J1023— 757 is equivalent to 19 pc at a distance of 6 kpc, which predicts 
a very fast mean expansion velocity of 4000 km s -1 over 4.6 kyr. The pulsar pseudo-luminosity 
distance places it much closer at 2.4 kpc (based on inferred beaming and gamma-ray efficiencies), 
though the scatter in inferred luminosities in radio-loud LAT pulsars translates to uncertainties in 
this estimate of the order of factors of 2-3 (|Saz Parkinson et al.ll2010l) . Both pulsar efficiency and 
PWN expansion velocity would be anomalously high at 6 kpc, so we adopt the pseudo-distance of 
2.4 kpc. At this distance the pulsar spin-down power (1. 1 x 10 37 erg s~ x ) can easily account for the 
VHE luminosity above 380 GeV of 1.4 x 10 34 d\ A erg s -1 . 



At longer wavelengths the vicinity of Westerlund 2 has undergone extensive study. Archival 
Chandra data indicate a faint source coincident with PSR J 1023— 5746, with an X-ray index of 
r = 1.2 ± 0.1 and unabsorbed 0.5 - 8 keV flux of l-3±o;| x 10~ 13 erg c m" 2 s" 1 , though th is does 
not affect modelling of the extended nebula. Recent Suzaku observations ( Fujita et aPl2009|) found 
no sign of diffuse non-thermal emission within the TeV contours, and placed a 0.7 — 2 keV upper 
limit on the diffuse flux from the entire XIS field of view of 2.6 x 10~ 12 erg cm -2 s _1 . Fujita 
and collaborators also note that it is unlikely that strong X-ray emission extends beyong this field 
since their upper limit is consistent with the one derive d using the wide HXD field (34' x 34'). 
Investigations of molecular clouds toward Westerlund 2 (IFukui et al]|2009i) show features of a few 
x 1O 4 M , though CO observations indicate a low density of gas (likely n < 1 cm -3 ) in the region 
that coincides with the bulk of the TeV emission. Radio observations of RCW 49 (the H II complex 
surro unding Westerlund 2) found a flux of 210 Jy at 843 MHz in the core (IWhiteoak & Uchida 
1997); this provides a non-constraining upper limit on the radio flux corresponding to the gamma- 
ray source. 

We computed SEDs (Spectral Energy Distributions) from evo lving electron popu lations over 
the lifetime of the pulsar in a series of time steps, as described in (|Abdo et alj|2010ch . As pulsars 
spin down, they dissipate rotational kinetic energy via 



e = inn 



with n the angular frequency and I the neutron star's moment of inertia, assumed to be 10 45 g cm 2 . 
This energy goes into a magnetized particle wind, and for magnetic dipole spindown of the pulsar 



n oc n 3 



(2) 
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Integrating equation 2 yields the age of the system ([Manchester & Taylorlll977|) : 

T 





(3) 



where P is the initial spin period, P the period derivative. For P <C P this equation reduces to 
the charactistic age of the pulsar r c = P/2P. The spin-down luminosity of the pulsar evolves as 
dPacini&Salvatilll973l) : 

E = E (1 + -)- 2 (4) 

T"0 



with the initial spin-down timescale defined as 



t 



2Pn 



(5) 



with P the initial spin period derivative. Given that the current P, P, and E are known, once 
an initial period is selected the age and spin-down history of the system is determined according 
to the equations above. We assume a particle dominated wind such that the wind magnetization 
parameter a ~ 1CT 3 . Therefore the power injected in the form of electron/positron pairs is E e = 
0.999E. 

As the distribution of particles expand with the PWN, they lose energy through adiabatic 
cooling, though synchrotron cooling typically dominates for the earliest phase of PWNe evolution. 
We assume that the radius R of the PWN scales linearly with time, and we select a magnetic 
field depende nce of B oc t~ L5 . B oth these behaviors closely mimic the behavior of B and R 



computed by iGelfand et al.1 (|2009|) for early stage PWN evolution prior to the compression and 
reexpansion phases caused by the interaction of the reverse shock. Selection of appropriate pho ton 
fields is crucial to accurate determination of IC fluxes. We therefore follow IPorter et al.l (I2006|) in 
estimating photon fields (CMBR, dust IR, and starlight) at the appropriate Galactic radii, unless 
local studies provide better estimates than these Galactic averages. 

To compute the PWN SED we inject at each time step a power-law spectrum of relativistic 
electrons with a high energy exponential cutoff. We also employ a low energy cutoff of 10 GeV 
for the electron spectrum , which is within the realm of minimum particle energies considered by 
Kennel & Coronitil (I1984Q . The energy content of this particle population varies with time follow- 
ing the pulsar spin down (eq. 4), though we treat the index and cutoff energies as static. We then 
adjust the size and magnetic field according to the models described above. Finally, we calculate 
the subsequent particle spectrum at time t + St by calculating the energy loss of the particles due 
to adiabatic losses as well as radiation losses from synchrotron and IC (including Klein-Nishina 
effects). Injection (and evolution) occurs in time steps much smaller than the assumed age. 
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Model fitting is achieve d by minim i zing the x 2 between model and data using the downhill 



simplex method described in Press et al.l (I1992() . We consider three variables: the initial spin pe- 
riod, electron slope, and high energy electron cutoff. With only an X-ray upper limit, the mean 
magnetic field within the gamma-ray source is poorly constrained, so we fix the current magnetic 
field to 5 fj,G (which is the best value obtained when we allow the magnetic field to vary), or ~ 2 
mG at pulsar birth. For each ensemble of these 3 variables we evolve the system over the pulsar 
lifetime and calculate x 2 ■ The simplex routine subsequently varies the parameters of interest to 
minizimize the fit statistic. We estimate parameter errors by computing x 2 for a sampling of points 
near the best fit values and using these points to fit the 3— dimensional ellipsoid describing the 
surface of A% 2 = 2.71. Under the assumption of Gaussian errors, the minima and maxima of this 
surface give the 90% errors of the parameters. 

For the assumed Galactic radius of PSR J 1023— 5746, dust IR photons typically peak at ~ 
T = 30 K with a dens ity ~ 1 eV cm~ 3 , while stellar photons peak at ~ T = 2500 K with a 



density 2 eVcm -3 ([Porter et al.ll2006|) . With these photon fields (and CMBR) we apply the 
model described above. Figure [3] indicates that IR photons dominate IC scattering above 10 GeV, 
with all three photon fields contributing for lower energies. For the best fit we find x 2 — 13.7 
for 8 degrees of freedom, with an electron power law index of 2.44 ± 0.06, high energy cutoff at 
60 ± 45 TeV, and initial spin period of 63 ± 17 ms. These parameters imply « 3 x 10 48 erg have 
been injected in the form of electrons, and an age of 3100 years. 



A hadronic origin for the observed gamma rays is also possible, and we follow IKelner et al 



(|20Q6Q in calculating the photons from proton-proton interactions and subsequent n° and r^-meson 
decay. Proton-proton interactions also yield 7r ± mesons which decay into secondary electrons, 
which we evolve in time. The t imescale for pion production via p-p interactions is given by r. 



pionpi 

1.5 x 10 8 (n/1 cm -3 ) -1 years (|Blumenthall 1 1 9701) : this timescale is significantly greater than the 



expected age of the system, so the proton spectrum is treated as static. We are able to fit the gamma- 
ray data only if the energy in protons exceeds 2 x 10 50 (n/1 cm" 3 ) d 2 .4 erg, with a x 2 ~ 15 for 
8 degrees of freedom. A hadronic origin for the gamma-rays is therefore energetically disfavored 
unless the gas density is much greater than 1 cm~ 3 throughout the bulk of the VHE emitting region. 
Yet we cannot rule out such an origin in the confused region around Westerlund 2, even though a 
PWN origin is reasonable given the fit parameters discussed above. 

Independent of the origin of the gamma rays, the lack of X-rays from the immediate vicinity of 
PSR J1023— 5746 is perplexing given its extremely high spin-down luminosity. One possibility is 
that electrons rapidly escape from the inner nebula into a low pressure bubble with correspondingly 
low magnetic field. For an electron conversion efficiency of ~ 1, at the current E after a mere rs 2 
years enough electrons are present in the inner nebula to recreate the observed X-ray flux for a 
20 /iG field appropriate for a termination shock. This timescale is comparable to the time for 
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particles to reach the termination shock. P ost-shock flow in PWNe, as determined by torus fitting, 
is typically w 0.7 c (iNg & Romanill2008|) : at this velocity particles will traverse the ~ 8" X-ray 
nebula surrounding J1023 in ~ 0.5^2.4 year. 



6. Discussion 



6.1. Constraints on pulsar modeling 



The high-quality statistics obtained with the Fermi-hAT both on the light curves and the spec- 
tra of the 54 pulsars detected allow a more detailed comparison with theoretical models than pre- 
viously possible. The detection or lack of significant emission in the off-pulse interval can also 
be used to discriminate between the different models. Currently, there are two classes of models 
that differ in the location of the emission region. The first compris es polar cap (PC) models whic h 
place the emission near the magnetic poles of the neutron star ( Daugherty and Hardin g[ 1996 ). 
The s econd class of outer magnetosphere models consists of the outer gap (OG) models (|Romani 
1996J), in which the emission e xtends between the null charge surface and the light cylinder, the 
two-pole caustic (TP C) models (IDvks an d Rudakl 120031) which might be realized in slot gap (SG) 



acceleration models (|Muslimov & Harding] 2004) . in which the emission takes place between the 
neutron star surface and the ligh t cylinder along the last open field lines, separatrix layer (SL) 
models (IBai & Spitkovskyl 120091) . in which emission takes place from the neutron star surface to 
outsid e the light cylinder and finally pair- starved polar cap (PSPC) models (|Muslimov & Harding 
20041) . where emission takes place throughout the entire open field region. Obs ervations by Fermi 
of simp le exponential cutoffs in the spectrum of Vela and other bright pulsars (lAbdo et al.ll2009al 
2010b|) . instead of super-exponential cutoffs expected in PC models, have clearly ruled out this 
class of model for Fermi pulsar emission. The outer magnetospher e models make different pre 



dictions for the level of of f -pulse emission. Classic OG models (|Romani & Yadigaroglulll995 



Chengl iRuderman & Zhangl : iRomani & Wattersll2010|) . for which there is no emission below the 
null charge surface, predict no off-pulse emission except at very small inclination angles and large 
viewing angles near 90 degrees. TPC models predict pulsed emission over m ost of the rotational 
phase at a level that dep ends on inclination, viewing angle, and gap width ( Venter et al.1 12OO9I : 
Romani & Watters 201o|) . In general, light curves for larger gap widths, expected for middle- 
aged and older pulsars in the SG model and when the viewing direction makes a large angle to 
the mag netic axis, have higher levels of off-pulse emission. The force-free magnetosphere SL 
model (IBai & Spitkovskyll2009|) also predicts light curves with off-pulse emission, since some ra- 
diation in this case also comes from below the null surface. PSPC m odels are expected to op erate 
in old and millisecond pulsars and predict off-pulse emission as well (iVenter & Harding Il2010|) . 
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Among the 54 pulsars analyzed in this paper, only 10 show a significant signal in their off- 
pulse, 7 of which are likely of magneto spheric origin (J0034-0534, J0633+1746, Jl 8 13- 1246, 
J1836+5925, J2021+4026, J2055+2539, J2124-3358). Two of the 7 showing off-pulse emission, 
J0034-0534 and J2124-3358, are millisecond pulsars. J1836+5925 with E = 1.2 x lO^ergs -1 , 



J2055+2539 with E = 5 x lO^ergs -1 and J0633+1746 (Geminga) with E = 3.3 x 10 34 



erg s 



-i 



have among the lowest spin-down luminosities of the normal Fermi detected pulsars. While 
J1813— 1246 and J2021+4026 have higher E (6.3 x lO^ergs" 1 and 1.1 x lO^ergs" 1 respec- 
tively), both have unusually wide gamma-ray pulses. 

As can be seen in the light curves presented in Figures |4] and [51 the level of off-pulse emission 
of these 7 pulsars greatly varies. The highest levels of off -pulse emission are found for J2021+4026 
with ~ 40%, J1836+5925 with ~ 35%, J0034-0534 and J2124-3358, with ~ 20% of the peak 
heights, while lower levels are found for J2055+2539 and J1813— 1246 with ~ 10% and Geminga 
with ~ 5% of the peak heights. In the case of TPC models, the highest off-pulse levels in light 
curves with two widely-spaced peaks ar e produced for inc lination angle a > 80° and viewing 
angle ( < 40° or a < 40° and ( > 80° J Venter et alJl2009h . In the case of OG models and wide 
peak separations, high le vels of off-pulse emission are produced only for a > 85° and £ < 30° 
(IRomani & Wattersll2010|) . For both types of model, a and £ must be very different (i.e., we are 
viewing the gamma-ray emission at a large angle to the magnetic axis) and these are precisely 
the conditions for which our line-of-sight does not cross the radio beam, and for which the pulsar 
should be radio quiet or radio-weak. In fact , all of the non-millisecond pulsars with si gnificant 
levels of off-pulse emission are radio quiet (|Saz Parkinson et al.l 12010c iRay et al.l l2010i) . In the 
case of the millisecond pulsars having large polar caps and small magnetospheres, the radio beams 
are thought to be much larger and a significant fraction of the gamma-ray beam size. Therefore, 
we may still view the radio beams at large angle from the magnetic pole. The light curve of 
J0034— 0534 shows two narrowly-spaced peaks which can be fit in both TPC and OG models 
for a = 30° and C = 70°, but off-pulse em ission is predicted in this case only for TPC models 
(Ubdo et alJboiOec Iventer & Harding I bold) . T he light curve of J 2124-3358 has actually been 
best fit with a PSPC model a = 40° and ( = 80° dventer et alJb009h . which also predicts off -pulse 
emission. In general, the detection of off-pulse emission in these pulsars constrains the outer gap 
solutions to a much greater degree than for TPC/SG or SL solutions. 



6.2. Constraints on pulsar wind nebulae candidates 

We searched for significant emission in the off-pulse window of 54 gamma-ray pulsars de- 
tected by Fermi-L AT and found only one convincing pulsar wind nebula candidate, J 1023-5746 
(besides the Crab Nebula and Vela-X). However, flux upper limits derived on the steady emis- 
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sion from the nebulae offer new constraints on sources already detected in the Te V range (e.g. the 
PWN e in the Kookaburra complex). Additionally, some PWNe were proposed by lBednarek & Bartosik 
J2005I) as promising sources of 7-ray emission in the GeV energy range, especially PSR J0205+6449 
and PSR J2229+61 14. We review some interesting cases in the following. 



6.2.1. PSR J0205+6449 and the PWN 3C 58 



The radio source 3C 5 8 was recognized early to be a supernova remnant (SNR G130.7+3.1) 
and classified as a PWN by IWeiler & Panagial (| 19781) . X-ray observations revealed a non-thermal 



spec trum with the photon index becoming steeper toward the outer region of the nebula (|Slane et al 
2004|). Flat spectrum radio emission S v oc v~ 0A2 covering r oughly 10' x 6' extends up to ~ 100 



GHz (Green 



(Slaneetal 



1986; Morsi& Reich 



20081) . and X-ray (ISlane et al. 



1987; 



Salte r et al.l ll989) and corresponds well with infrared 



2004) morphologies. Subsequent Chandra X-ray Ob- 



servatory observations detected the central pulsar of 3C 58, PSR J0205+6449. The pulsar has a 
very high spin-down power of 2.7 x 10 37 erg s" 1 and a characteristic age of 5400 yr. 3C 58 has 
often been associated with S N 1181 (IStephenson et al.ll2002|) . However, recent investigations of 
the dynamics of the system (iChevalierl |2005|) . and the velocities of both the radio expansion and 
optical knots imply an age of ~ 2500 yr, closer to the characteristic age of PSR J0205+6449. 
At TeV energies, both the VERITAS and MAGIC telescopes o bserved this source and did not fin d 
any evidence for 7-ray emission at the position of the pulsar (|Anderhub et al.ll201Cc lAliul 12008). 
The upper limits derived from their observations are consistent with the Fermi upper limits ob- 
tained in the 100 MeV - 100 GeV energy range of < 12.9 x 10~ 12 erg cm~ 2 s _1 . This upper limit 
implies a non-constraining 100 MeV - 100 GeV efficiency of < 4 x 10~ 4 - 6 x 10~ 4 at a distance 
of 2.6-3.2kpc. 



6.2.2. PSR J0633+1746 - Geminga 



The Geminga pulsar is the first representative of a population of radio-quiet gamma-ray pul- 
sars, and has been in tensely studied since its discovery a s a ga mma-ray source by SAS-2, more 
than thirty years ago (lFichtel et al.l (|1975 ); Knif fen et al.l d 1975k The subsequent ROSAT detec- 
tion of periodic X-rays from this source (IHalpern & Holtlll992l) prompte d a successful search for 
periodicity in high-energy gamma-rays with EGRET (IBertsch et al.ll 19921) X-ray observations with 
XMM-Newton and Chand ra observations ind icate a highly structured pulsar wind nebula extend- 
ing ~ 50' from the pulsar (|Pavlov et al.ll2010|) . 



The Milagro Collaboration recently reported a 3.5a source coincident with Geminga of extent 
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~ 2.6°, likely the result of a PWN (lAbdo et al.ll2009er) . While VHE emission from the vicinity 
of this unique pulsar is certainly of interest, the Fermi upper limits on off-pulse emission above 
10 GeV are not very constraining to SED modeling of the single VHE point. The lack of data, 
compounded by the lack of X-ray data covering the scale of the Milagro source, renders SED 
modeling uninformative. 



6.2.3. PWNe in the Kookaburra complex: the Rabbit and K3 



The complex of compact and extended radio sources called Kookaburra (Ro berts et al.lll999|) 
covers nearly one square degree along the Galactic plane around I = 313.4°. This region has 
been extensively studied to understand the nature of the unidentifie d EGRET source 3EG J 1420- 
6038 (IHartmanetal.lll999r) . In the North wing of the Kookaburra. Id' Amico et all (1200 lb discov- 
ered the radio pulsar PSR J 1420— 6048, a young energetic pulsar with period 68 ms, characteristic 
age r c = 13 kyr and spin-down power E = 1.0 x 10 37 erg s _1 . Subsequent ASCA and Chan- 
dra observations revealed an X-ray nebula surrounding the pulsar called K3. In the South- Western 
wing, a bright X-ray emission called the Rabbit has been proposed as a plausible PWN contributing 
to the 7-ray signal detected by EGRET. At TeV energies, the H.E.S.S. co llaboration reported th e 
detection of two bright sources coincident with the Kookaburra complex (|Aharonian et al.ll2007|) . 
HESS J1420-607 is centered just North of PSR J1420-6048, with best fit position overlapping 
the pulsar position. HESS J1418— 609 appears to correspond to the Rabbit nebula. However, the 
gaps in spectral coverage between H.E.S.S. and EGRET did not allow a clear statement if the 
EGRET source is really associated with the X-ray and TeV emission. Finally, the 7-ray detec- 
tion by Fermi-LAI of the radio loud pulsar PSR J1420— 6048 and the discovery of a radio-quiet 
PSR J1418— 6058, likely powering the Rabbit PWN, brought a new light on this region. 
In this paper, we searched for 7-ray emission in the off-pulse of both pulsars, PSR J1420— 6048 
and PSR J141 8—6058, and did not detect any significant sign al. The 95% CL upper l imits reported 
in Tables [2] and [3] are consistent with the prediction made by IVan Etten & Romanil (|2010 ) that do 
not expect any detection by Fermi before 10 yr of observations in the most optimistic scenario for 
PSR J1420— 6048. The VHE spectrum of the Rabbit nebula is very similar to that of K3, so sim- 
ilar emission mechanisms from the Rabbit would correspondingly predict little chance of Fermi 
detection. 



6.2.4. PSR J1833-1034 and G21.5-0.9 



G21.5-0.9 was classified as one of abou t ten Crab-like SNR and predicted in 1995 to be a 
gamma-ray source dde Jager & Hardingll 19951) . Chandra observations revealed the composite na- 
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ture of the remnant, c onsisting of a centrally peaked PWN and a 4' shell (|Bocchino et all 12005 



Safi-Harbet all 1200 II) . The 61.8 ms pulsar PS R J1833-10 34 powering the PWN was discov- 



ered recently through its faint radio emission (ICamilo et al.ll2006() . With a spin-down power of 
E = 3.3 x 10 37 erg s _1 , PSR Jl 833- 1034 is one of the most energetic pulsars in the Galaxy. INTE- 



GRA L observations revealed that the PWN is also bright in the hard X-ray regime (Ide Rosa et al. 



2009|). At TeV energies, G2 1.5-0. 9 was rece ntly detected by H.E.S.S . and has a point-like na- 



ture and a hard spectrum (T = 2.08 ± 0.22) (|Djannati-Atai et al.ll2007|) . The flux of this source 
is only 2% of that of the Crab Nebula. Although the spin-down age of this extremely energetic 
pulsar is 4.6 kyr, VLA me asurements of the PWN expansion speed place the age at 870i°° 50 yr 
( Bietenholz & Bartel 12008). Kinema tic H I and CO distance measurements place G21.5-0.9 some 
4.7 kpc distan t jTian & Leahy 20081) . The lack of detection is consistent with previous modeling 
undertaken by Ide Jager et al.1 (I2008|) . which predicted GeV gamma-ray flux well below the Fermi 
upper limits. 



6.2.5. PSR J 1907+0602 and its TeV PWN 



The TeV source M GRQ J 1908+06 was discovered by MILAGRO at a median energy of 20 
TeV (lAbdo et al]|2007|) with a flux ~80% of the Cr ab at these energies. It was subsequentl y de 



tecte d in the 300 GeV-20 TeV range by the H.E.S.S. (lAharonian et alJl2009h and VERI TAS ([Ward 



2008) experiments. The Fermi discovery of the radio-quiet pulsar PSR J1907+0602 (|Abdo et al 



2009d|) within the extent of the TeV source strongly suggests that HESS J 1908+063 is the PWN of 
PSR J1907+0602. 

The upper limits derived on i ts off-pulse emission (see Table |2]) are consistent with those previ- 
ously reported in lAbdo et al.l (|2010gb assuming a source extension of 0.3° and strongly suggest 
that the spectrum of HESS J 1908+063 has a low-energy turnover between 20 GeV and 300 GeV. 
The pulsar is offset from the H.E.S.S. centroid by 15', and assuming this marks the pul sar birthsite 
implies a velocity of ~ 400kms _1 for a distance of 3.2 kpc and an age of 20 kyr. lAbdo et al 



estimated an upper limit on the 2-10 keV X-ray flux from the pulsar and any arcminute- 
scale nebula of 1CT 13 erg cm -2 s _1 , though given the H.E.S.S. extension is 20' we cannot constrain 
the magnetic field within the larger nebula. Therefore without any data outside the VHE regime 
we can only conclude that a cooling break appears around 1 TeV. Determining whether that break 
occurs from old electrons in a low magnetic field or younger electrons in a high magnetic field 
requires either an X-ray detection or upper limit on the extended region. 
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6.2.6. PSR J2032+4127 and TeV 2032+4130 



TeV J2032+4130 was the first unclassi fied TeV source, initially detected by HEGRA (lAharonian et al 



20021) . and later confirmed by MAGIC (1 Albert et all 120081) . Subsequent XMM observations re- 
vealed a faint diffuse X-ray stru cture centered on the position of TeV J2032 +4130 with the same 
6' extension as the VHE source (IHorns et al.ll2007|) . Within the TeV error box lAbdo et alj (I20 09d) 



disco vered PSR J2032+4127, later confirmed by the radio detection of the pulsar (ICamilo et al. 



20091), hinting at a PWN origin for the TeV emission. The limited radio data and large errors on 
the XMM X-ray and VHE data prevent precise modeling of all PWN parameters, however, and the 
Fermi upper limits are approximately an order of magnitude too high to be constraining. 



6.2.7. PSR J2229+6114 and the Boomerang PWN 



The pulsar PSR J2229+6114 is as young as the Vela pulsar (characteristic age r c = 10 kyr), 
as energetic (E = 2.2 x 10 37 erg s _1 ) and powers the small boomerang- shaped radio and X - 
ray emitting PWN which is part of the supernova remnant G1 06.3+2.7 (IHalpern et al.ll2001b|) . 
The pulsar distance estimated from X-ray absorption is ~3 kpc (Halper n et al]|2001bl). w hile the 
dispersion measure us ed in conjunction with the NE2001 model ( Cordes & Laziol 2002|) yields 
a distance of 7.5 kpc ( AbdoetaL 2009cl): polarization s tudies and velocity maps of HI and CO 
emission imply a distance of 800 pc (|Kothes et al. 2006 ). In this paper, we have used the range 
0.8-6.5 kpc as reported in Table [T] 

At TeV energies, the MAGIC collaboration place d a constraining point- source upper limit of 



3.95 x 10 12 cm 2 s 1 at the position of the pulsar (lAnderhub et al.ll2010|) . However, recently. 



both the MILAGRO and the VERITAS collaboration reported a significant detection in this region. 
The centroid of the extended emission (full angular exte nt of 0.6° by 0-4°) d etected by VERITAS 
is located 0.4° away from the pulsar PSR J2229+6114 (lAcciari et al]|2009|) . The signal detected 
by MILAGRO is spread over a broad ~ 1° area encompassing the pulsar position and the main 
bulk of the remnant, which does not allo w a definitive association with a particular region of the 
SNR/pulsar complex (Ubdo et al.ll2009eh . 



Using 16 months of Fermi-LAI data, we do not report any significant detection at the position 
of the pulsar using its off-pulse photons. The upper limits deri ved are compatible wi th the non- 
detection reported by the MAGIC collaboration at TeV energies ( Anderhub et al.ll2010|) . An offset 
between the pulsar J2229+6114 and its PWN is thus required if the emission detected at TeV 
energies by VERITAS and MILAGRO is produced by the pulsar wind nebula. This is the case 
for several PWNe already detected by Cherenkov telescopes, such as HESS J 1825— 137 and Vela- 
X, and can be explained by the supernova explosion that occured in an inhomogeneous medium 
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leading to an asymmetric reverse shock that displaced the PWN (|Blondinll2001|) towards lower 
densities. However, the TeV emission detected by VERITAS is coincident with the location of 
molecular clouds which disfavors such a scenario and supports a hadronic origin where 7-rays are 
produced via proton-proton interactions within the molecular clouds. 



6.3. Population study of pulsar wind nebulae as observed by Fermi-LAT 

Among the large sample of pulsars analyzed in this paper, two sources are firmly identified 
as PWNe (the Crab Nebula and Vela-X), one source is suggested as a highly plausible candi- 
date (the emi ssion in the off-puls e of PSR J1023— 5746), in addition to the pulsar wind nebula in 



MSH 15— 52 lAbdo et al.1 (|2010h|) . As can be seen in Figure [6l the pulsars powering these PWNe 
are all young (in the range 1 to 10 kyr) and bright (E > 7 x 10 36 ergs" 1 ); their associated PWNe 
are detected by Cherenkov telescopes in the TeV range. Remarkably, these 4 PWNe candidates 
have a low 7-ray efficiency in the Fermi-LAT energy range, Lpi Y N < 1CT 2 , Vela-X, Crab and 
MSH 15—52 having even lower Lp W N (see Figure [7]). This implies that most of the 44 upper 



E 

limits derived using 16 months of LAT observations are not yet constraining. However for three 
objects (J0659+1414, J1833-1034 and J0205+6449) the upper limits in Figure Hare well below 
2 x lCT 3 -^ suggesting that some pulsars are less efficient at producing GeV PWN flux than J1023- 
5746 and B 1509-5 8. Table [2] indicates three potential pulsar wind nebulae candidates at 4<r level: 
J0007+7303, J1028-5819, J1709-4429. More data are needed to confirm the detection these 3 
sources and significantly probe the variation of the PWN flux with the spin-down power or with 
the ambient medium. One can also note two pulsars with significant emission in their off-pulse 
window and very high efficency: J1836+5925 and J2021+4026. These two pulsars show a simi- 
lar gamma-ray luminosity in the off-pulse and in the on-pulse intervals (renormalized to the total 
phase interval); their luminosities are greater than the spin-down power, assuming a uniform phase- 
averaged beaming across the sky. However, the association distances for the gamma-ray selected 
pulsars must be treated with caution. 



Reporting the Fermi detection of HESS J1640— 465. 1Slane et al.l (|2010|) invoked an excess of 
low energy electrons to account for the GeV data and to explain the continuity in photon spectral 
index between the Fermi and HESS energy bands. While many of the upper limits on PWNe flux 
catalogued here are insufficiently constraining to rule out such a situation, we also see no evidence 
in support of such an electron spectrum for any object in t he catalog. We also see no indication of 



dual electron populations, as proposed to explain Vela-X (lAbdo et al.H2010q) . The GeV and TeV 
peaks in Vela-X are comparable in flux, and this would put the GeV flux very near the Fermi sensi- 
tivity for a number of Vela-like pulsars with TeV counterparts: PSRs J1418-6058, J1420-6048, 
J1709— 4429, and J1907+0602. Indeed, if most PWNe boast simple electron populations with a 
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power-law index of ~ 2.3, the majority of these PWNe may prove elusive to Fermi. 
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Table 1. Observatories, off -pulse definitions and distances of the 54 pulsars analysed 



PSR 



ObsID Off-pulse definition Distance (kpc) Observation period rejected (MJD) 



J0007+7303 9 


L 


0.4 - 0.8 


1.4±0.3 


J0030+O451 


N 


0.7-1.1 


0.300±0.090 


J0034-0534 


N 


0.45 - 0.85 




J0205+6449 9 


G,J 


0.7- 1.0 


2.6-3.2 


J02 18+4232 


N 


0.9- 1.1 


2.5^1 


J0248+6021 9 


N 


0.7- 1.1 


2.0±0.2 


J0357+32 


L 


0.35 - 0.85 




J0437-4715 


P 


0.7- 1.2 


0.1563±0.0013 


J0534+2200 


N, J 


0.5 - 0.85 


2.0±0.5 


J0613-0200 


N 


0.6- 1.05 


0.48t°;J? 


J0631+1036 9 


N, J 


0.9- 1.15 


0.75-3.62 


J0633+0632 


L 


0.6-0.8 




J0633+1746 


L 


0.67 - 0.87 


U.ZJU_g 062 


J0659+1414 


N, J 


0.45 - 1.0 


n 9RC+0.033 
U.ZSS_ 27 


J0742-2822 9 


N, J 


0.8-1.4 


1 fy7+l-38 
z - u/ -1.07 


J0751+1807 


N 


0.7- 1.05 




J0835-4510 


P 


0.7- 1.0 


0.2871°;°!? 


J1023-5746 9 


L 


0.85- 1.13 


2.4 


J1028-5819 


P 


0.8-1.05 


2.33±0.70 


J1044-5737 


L 


0.75-1.1 


1.5 


J1048-5832 


P 


0.7-1.05 


2.71±0.81 


J1057-5226 


P 


0.7-0.2 


0.72 ±0.2 


Jl 124-5916 9 


L 


0.92 - 0.08 


48 +0.7 
t.O_ 1 2 


J1413-6205 9 


L 


0.7-0.15 


1.4 


J1418-6058 


L 


0.55 - 0.90 


2-5 


J1420-6048 


P 


0.6- 1.1 


5.6±1.7 


J1429— 5911 


L 


0.85-0.1 


1.6 


J1459-60 


L 


0.34 - 0.69 




J1509-5850 


P 


0.6-1.0 


2.6±0.8 


J1614-2230 


G 


0.92- 1.14 


1.27±0.39 


J 1709 -4429 9 


P 


0.65-1.1 


1.4-3.6 


J1718-3825 


N, P 


0.65- 1.15 


3.82±1.15 


J1732-31 


L 


0.54 - 0.89 




J1741-2054 


L 


0.67- 1.18 


0.38±0.11 


J1744-1134 


N 


0.15-0.35 


A i<; 7 +0.043 


J1809-2332 


L 


0.45 - 0.85 


1.7±1.0 


J1813-1246 9 


L 


0.72 - 0.84 




J1826-1256 


L 


0.60 - 0.90 




J1833-1034 


G 


0.75-1.1 


4.7±0.4 


Jl 836+5925 


L 


0.16-0.28 


<0.8 


J1846+0919 


L 


0.65 - 1.0 


1.2 


J1907+06 


L 


0.51-0.91 




J1952+3252 


J, N 


0.7- 1.0 


2.0±0.5 


J1954+2836 


L 


0.85 - 0.2 


1.7 


J1957+5033 


L 


0.6 - 0.05 


0.9 


J1958+2846 


L 


0.55 - 0.90 





54870 - 54940 
55161 -55181 



54682 - 54743 



55084-55181 
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Table 1 — Continued 



PSR 


ObsID 


Off-pulse 


definition 


Distance (kpc) Observation period rejected (MJD) 


J2021+3651 9 


G 


0.75 


- 1.05 


2 1+ 2 ' 1 
z - 1 -1.0 


J2021+4026 


L 


0.16 


- 0.36 


1.5±0.45 


J2032+4127 


L 


0.30 - 0.45 & 0.90 - 0.05 


1.6-3.6 


J2043+2740 


N, J 


0.68 


-0.08 


1.80±0.54 


J2055+2539 


L 


0.6 


-0.1 


0.4 


J2124-3358 


N 


0.1 


-0.5 


25+ - 25 


J2229+6114 9 


G, J 


0.68 


- 1.08 


0.8-6.5 


J2238+59 


L 


0.65 


- 0.90 





Note. — Column 1 lists the pulsars; a "g" indicates that one or several glitches occured during the observation 
period. For some pulsars, these glitches led to restrict the dataset to avoid any contamination of pulsed emission 
during the glitch: the observation period rejected in these cases is indicated in column 5 (Modified Julian Day). 
Column 2 indicates the observatories that provided ephemerides: "G" - Green Bank Telescope; "J" - Lovell telescope 
at Jodrell Bank; "L" - Large Area Telescope; "N" - Nancay Radio Telescope; "P" - Parkes radio telescope. Column 
3 lists the off-pulse phase range used in the spectral analysis. Column 4 presents the best known distances of 54 the 
pulsars analyzed in this paper. 



Table 2. Spectral fit results for 54 LAT-detected pulsars 



PSR 


TS 


^0.1 — 100 
(10~ 9 ph cnr 2 s _1 ) 


^0. 1 — 100 
(10~ 12 erg cm~ 2 s _1 ) 


r 


T uminncitv 

HO 33 ere-s -1 ) 


J0007+7303 


24.3 


<63.23 


<69.94 




< 16.40 


J0030+0451 


3.4 


<7.07 


<7.83 




<0.08 


J0034-0534 


29.1 


17.26 ± 5.70 


11.09 ± 2.68 


2.27 ±0.17 


25+ ' 75 
—0.25 


J0205+6449 


1.3 


< 11.63 


<12.88 




< 10.42- 15.78 


J02 18+4232 


1.2 


<12.33 


<13.65 




<10.21 -26.13 


J0248+6021 


0.3 


<7.77 


<8.59 




<4.11 


J0357+32 


0.0 


<3.94 


<4.36 






J0437-4715 


10.7 


<8.50 


<9.41 




<0.03 


J0534+2200 a 


2775.6 


980.00 ± 70.00 


540.92 ± 46.73 


2.15 ±0.03 


258.88±151.81 


J0613-0200 


4.0 


<6.74 


<7.46 




<0.21 


J0631+1036 


2.5 


<18.72 


<20.72 




< 1.39 -32.49 


J0633+0632 


6.3 


<32.50 


<35.97 






J0633+1746 


5101.2 


1115.54 ± 32.31 


749.44 ± 22.24 


2.24 ± 0.02 




J0659+1414 


0.6 


<5.04 


<5.58 




<0.05 


J0742-2822 


0.0 


<5.99 


<6.63 




<3.40 


J0751+1807 


6.4 


<9.52 


<10.53 




<0.45 


J0835-4510 6 


284.3 


405.44 ± 26.75 


210.25 ± 13.87 


2.30 ±0.10 


2 07+ ' 41 

— 0.38 


J1023-5746 


25.1 


1.33 ± 1.14 


27.58 ± 13.73 


1.05 ±0.36 


19.01 ± 9.46 


J1028-5819 


20.8 


<88.79 


<98.27 




<63.83 


J1044-5737 


0.0 


< 11.93 


<13.20 




<46.05 


J1048-5832 


0.0 


<15.33 


< 16.96 




< 14.90 


J1057-5226 


1.2 


< 10.26 


<11.35 




<0.70 


J1124— 5916 


0.0 


< 12.09 


<13.38 




<36.88 


J1413-6205 


5.3 


<4.83 


<5.34 




< 14.72 


J1418-6058 


10.3 


<77.73 


<86.03 




<41.17- 257.33 


J1420-6048 


3.0 


<125.98 


< 139.42 




<523.13 


J1429— 5911 


0.0 


<21.26 


<23.53 




<88.29 


J1459-60 


1.2 


<23.17 


<25.64 






J1509-5850 


1.1 


<25.47 


<28.19 




<22.80 


J1614-2230 


5.3 


<22.01 


<24.36 




<9.55 


J1709-4429 


16.5 


<35.59 


<39.39 




<61.08 


J1718-3825 


0.0 


<8.53 


<9.44 




< 16.48 


J1732-31 


0.0 


<7.40 


<8.19 






J1741-2054 


0.3 


< 10.52 


<11.64 




<0.20 



Table 2 — Continued 



PSR 


TS 


Fo. l-ioo 


Go. i-ioo 


r 


Luminosity 






SI r\ Q 1 9 1 \ 

(10 9 phcm 2 s 1 ) 


(10 -12 erg cm~ 2 s _1 ) 




(10 33 ergs -1 ) 


T1 744—1 114 


6.9 


<27 68 


<30 64 




<0 47 


Tl RfiQ— 9319 


1 9 


< 19 19 


<21 25 




<7 35 


T1R11 194ft 


38 1 




1 1 Q fla _l_ n 90 


Z.UJ 31 \). 1 + 






9 7 


<145 17 


<160 67 






T 1 1 014 





<9 38 


<10 38 




<27.43 


J 1 OJOTJ7ZJ 


99cn 6 


S7Q An -t- 98 S6 


S49 1 6 4- 34 fH 


9 07 -t- 03 
z.u / in u.uj 


9/; 77-1-1 91 


Tl S4A+091 9 


0.0 


<4.79 


<5.30 




< 10.66 




0.7 


< 17.02 


< 1 8.83 






J1952+3252 


2.4 


<16.88 


<18.68 




<8.94 


J 1954+2836 


2.4 


<21.49 


<23.78 




<99.04 


J1957+5033 


0.3 


<5.45 


<6.04 




<7.40 


J1958+2846 


2.6 


< 15.43 


< 17.07 






J202 1+3651 


15.8 


<91.48 


< 101.24 




<53.42 


J202 1+4026 


2229.1 


1603.0 ± 11.2 


888.12 ±8.56 


2.36 ± 0.02 


198.45±119.83 


J2032+4127 


1.2 


< 154.91 


<171.45 




<52.51 -265.86 


J2043+2740 


0.0 


<2.71 


<2.99 




< 1.16 


J2055+2539 


36.7 


38.41 ± 10.10 


17.59 ± 3.34 


2.51 ±0.15 


2.87 ±1.44 


J2124-3358 


64.6 


22.78 ± 6.43 


21.81 ±4.44 


2.06 ±0.14 


10+ ' 22 


J2229+6114 


0.4 


<14.05 


<15.55 




<1. 19-78.61 


J2238+59 


0.0 


< 14.92 


<165.10 







"The spectral parameters of the Crab Nebula are derived using Ubdo et alj feoiOot) . 

"The spectral parameters are derived assuming a uniform disk morphology as described in lAbdoetaljfcoiOch . 



Note. — Results of the maximum likelihood spectral fits for the off-pulse emission of LAT gamma-ray pulsars (see 
Section[4} between 100 MeV and 100 GeV. Pulsar wind nebula spectra are fitted with a power-law model (photon index 
r, photon flux F and energy flux G) assuming a point-source at the position of the pulsar. The test statistic (TS) for 
the source significance is provided in Column 2, the photon flux F and the energy flux are reported in Columns 3 and 
4 while the photon index is listed in Column 5 when TS > 25. The photon flux and energy flux obtained from the 
likelihood analysis are replaced by a 2a upper limit when TS < 25 (assuming a photon index T=2). The total gamma- 
ray luminosity L 7 is listed in column 6. The error on the photon flux and the photon index only include statistical 
uncertainties while the error on include the statistical uncertainties on the flux and the distance uncertainties. 



Table 3. Spectral fit results for 54 LAT-detected pulsars 



HK 




n 1 1 rk»v 

U. 1 - 1 Uc V 






1 - 1U Uc V 






1 n 1 on n^v 

1 U - 1 UU Uc V 






l o 


^0.1—1 

pn cm s ) 


p 




^1 — 10 
(_iu pn cm s ) 


p 


l o 


-FlO— 100 
{±\J pn cm s ) 


p 


JO0O7+7303 


22.8 


<54.39 




9.2 


<2.75 




0.0 


<0.24 




J0030+0451 


6.3 


<17.11 




3.7 


<0.67 




0.0 


<0.18 




J0034-0534 


26.9 


12.36 ±7.15 


1.47 ±0.60 


16.5 


0.97 ± 0.35 


3.14 ±0.78 


0.0 


<0.17 




J0205+6449 


0.0 


< 12.35 




2.5 


<1.99 




1.6 


<0.30 




J02 18+4232 


0.3 


< 24.46 




1.1 


<1.39 




0.0 


<0.47 




J0248+6021 


0.0 


<5.40 




0.0 


<0.35 




0.3 


<0.25 




J0357+32 


().() 


<9.15 




0.0 


<0.56 




0.0 


<0.14 




J0437-4715 


4.1 


<13.83 




8.7 


<0.94 




0.0 


<0.18 




J0534+2200 a 


1054.5 


785.14 ±45.37 


3.20 ± 0.07 


1206.9 


22.93 ± 1.44 


1.59 ±0.10 


830.7 


5.12 ±0.56 


1.91 ±0.19 


J0613-0200 


5.5 


<26.97 




0.5 


<0.72 




0.0 


<0.16 




J0631+1036 


4.6 


<49.87 




0.27 


<0.18 




0.0 


<1.11 




J0633+0632 


7.4 


<67.47 




1.3 


<2.85 




0.0 


<0.42 




J0633+1746 


3377.3 


837.66 ± 32.20 


1.81 ±0.05 


2028.4 


65.41 ± 3.08 


3.26 ±0.11 


0.0 


<0.35 




J0659+1414 


0.5 


< 11.94 




1.7 


<0.59 




0.0 


<0.13 




J0742-2822 


0.0 


<24.54 




0.0 


<0.94 




0.0 


<0.13 




J0751+1807 


1.7 


< 13.99 




11.4 


<1.46 




0.0 


<0.19 




J0835-4510 6 


199.8 


329.76 ± 34.54 


2.15 ±0.11 


97.9 


18.36 ± 2.33 


2.22 ± 0.20 


2.4 


<0.89 




J1023-5746 


0.0 


<12.55 




0.9 


<2.42 




17.2 


0.46 ± 0.22 


1.02 ±0.73 


J1028-5819 


15.8 


< 180.07 




15.0 


<7.18 




0.0 


<0.53 




J1044-5737 


0.6 


<45.53 




0.0 


<1.46 




0.0 


<0.36 




J1048-5832 


0.0 


<40.52 




0.4 


<2.12 




0.0 


<0.33 




J1057-5226 


0.8 


<22.44 




2.4 


<1.33 




0.0 


<0.21 




Jl 124-5916 


0.0 


<49.42 




0.2 


<2.52 




0.0 


<0.49 




J1413-6205 


2.8 


<55.63 




11.8 


<6.77 




0.0 


<0.37 




J1418-6058 


5.2 


<94.24 




6.9 


< 10.28 




1.1 


<0.66 




J1420-6048 


6.7 


< 29 1.45 




0.0 


<9.62 




1.1 


<0.83 




J1429— 5911 


0.0 


<54.23 




0.0 


<3.25 




0.0 


<0.48 




J1459-60 


5.0 


<68.17 




0.1 


<1.76 




1.5 


<0.56 




J1509-5850 


0.7 


<65.47 




0.4 


<3.33 




0.0 


<0.33 




J1614-2230 


2.5 


<34.37 




7.7 


<2.30 




0.0 


<0.43 




J1709-4429 


15.5 


<96.67 




3.3 


<2.36 




0.0 


<0.22 




J1718-3825 


0.5 


<56.84 




0.0 


<0.92 




0.1 


<0.24 




J1732-31 


0.0 


<35.79 




0.0 


<1.09 




0.0 


<0.25 





Table 3 — Continued 



PSR 




0.1 - 1 GeV 






1 - 10 GeV 






10 - 100 GeV 




TS 


-Fo.l-l 
(10 -9 ph cm -2 s -1 ) 


r 


TS 


-Fi-io 
(10" 9 ph cm- 2 s- 1 ) 


r 


TS 


Fio— ioo r 
(10~ 9 ph cm- 2 s- 1 ) 


J1741-2054 


0.9 


<27.77 




4.4 


<1.40 




0.0 


<0.14 


J1744-1134 


10.3 


<71.71 




4.9 


<1.92 




0.0 


<0.45 


J1809-2332 


2.8 


<50.50 




8.8 


<2.40 




1.2 


<0.27 


J1813— 1246 


32.7 


261.21 ± 73.15 


2.25 ± 0.27 


16.3 


8.43 ± 2.57 


3.05 ± 0.69 


3.6 


<1.09 


J1826-1256 


9.5 


<251.30 




3.4 


<5.81 




0.1 


<0.38 


J1833— 1034 


0.0 


<13.33 




0.1 


<1.67 




0.2 


<0.40 


J 1836+5925 


1381.5 


401.84 ±27.39 


1.56 ±0.09 


1014.1 


51.36 ± 3.89 


2.93 ±0.16 


0.0 


<0.74 


J 1846+09 19 


0.0 


<17.79 




0.0 


<0.61 




0.0 


<0.19 


J 1907+06 


1.9 


<70.95 




0.9 


<2.33 




0.0 


<0.20 


J 1952+3252 


1 .4 


<50.81 




1.3 


< 1 .58 




0.0 


<0.23 


J1954+2836 


1.4 


<47.98 




2.6 


<2.55 




0.0 


<0.25 


J1957+5033 


1.2 


<16.35 




0.2 


<0.57 




0.0 


<0.17 


J1958+2846 


0.0 


<27.14 




5.3 


<1.78 




0.0 


<0.26 


J2021+3651 


17.7 


85.90 ± 30.02 


1.90 ±0.31 


13.2 


<4.50 




0.0 


<0.25 


J202 1+4026 


1718.2 


1344.75 ±55.56 


2.03 ± 0.05 


936.2 


73.76 ± 3.93 


3.04 ±0.11 


12.16 


<1.24 


J2032+4127 


3.5 


< 133.39 




0.0 


<2.08 




1.3 


<0.56 


J2043+2740 


0.0 


<9.73 




0.0 


<0.76 




0.0 


<0.17 


J2055+2539 


35.3 


16.06 ± 10.90 


1.23 ±0.76 


23.3 


1.53 ±0.42 


4.89 ± 0.75 


0.0 


<0.13 


J2124-3358 


16.0 


16.75 ± 12.17 


1.83 ±0.70 


56.6 


2.41 ± 0.54 


2.34 ± 0.37 


0.0 


<0.21 


J2229+6114 


4.2 


<49.68 




0.0 


<1.36 




0.0 


<0.28 


J2238+59 


2.5 


<55.91 




0.0 


<1.54 




0.0 


<0.38 



a The spectral parameters of the Crab Nebula are derived using lAbdoetaljf2010dh . 

The spectral parameters are derived assuming a uniform disk morphology as described in lAbdo et all feoiOch . 



Note. — Results of the maximum likelihood spectral fits for the off-pulse emission of LAT gamma-ray pulsars (see Section[4}. The off-pulse spectra were fit with 
a power-law model (photon index V and photon flux F) assuming a point-source at the position of the pulsar. The results for the fits in the three energy bands are 
reported. The test statistic (TS) for the source significance is provided in Columns 2 (0.1 - 1 GeV)), 5 (1 - 10 GeV) and 8 (10 - 100 GeV). The photon flux F for 
each energy band is reported in Columns 3, 6, and 9; it is replaced by a 2<r upper limit when TS < 25 (assuming a photon index T=2). Columns 4, 7 and 10 list the 
photon index T for each energy band when TS > 25. Only statistical uncertainties are reported on the photon flux and the photon index. 



Table 4. Spectral fitting of pulsar wind nebula candidates with low energy component. 



PSR 


Go. i-ioo 


r 


-^cutoff 


TS cu t G 




(10~ 12 erg cm -2 s _1 ) 




(GeV) 




J0034-0534 


7.33 ± 2.01 ± 1.30 


0.62 ± 1.05 ±0.27 


0.7 ±0.48 ±0.10 


9.0 


J0633+1746 


544.01 ± 13.91 ± 54.58 


1.51 ±0.06 ±0.12 


1.41 ±0.14 ±0.09 


247.2 


J1813-1246 


116.24 ± 22.92 ±79.28 


2.65 ±0.14 ±0.26 




1.2 


Jl 836+5925 


349.64 ± 16.04 ±28.05 


1.33 ±0.10 ±0.06 


1.60 ±0.25 ±0.04 


99.8 


J202 1+4026 


737.14 ± 21.77 ± 125.06 


1.87 ±0.06 ±0.20 


2.24 ±0.37 ±0.51 


110.2 


J2055+2539 


12.23 ±6.14 ±6.09 


0.30 ± 1.40 ±0.69 


0.43 ±0.31 ±0.07 


22.4 


J2124-3358 


13.27 ± 3.02 ± 2.77 


0.88 ± 0.74 ± 0.34 


1.71 ± 1.06 ± 0.59 


10.4 



Note. — Results of the maximum likelihood spectral fits for pulsars showing a significant signal 
in their off-pulse at low energy. The fits used an exponentially cutoff power-law model with the 
energy flux Go.i-ioo> photon index T and cutoff energy E CVLto f[ given in columns 2, 3 and 4. The first 
errors represent the statistical error on the fit parameters, while the second ones are the systematic 
uncertainties as discussed in section I5T1 The significance of an exponential cutoff (as compared to a 
simple power-law) is indicated by TS cut oS in column 5. A value TS CVLto s < 9 indicates that the two 
models are comparable and we report the fit parameters assuming a simple power-law model. 
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Fig. 1. — Spectral energy distributions of the off -pulse emission of J0034— 0534 (top left), 
J0633+1746 (top right), J1813— 1246 (bottom left) and J1836+5925 (bottom right), renormalized 
to the total phase interval. The LAT spectral points are obtained using the maximum likelihood 
method described in section I5T1 into 6 logarithmically- spaced energy bins. The dotted-dashed 
green line presents the result obtained by fitting a power-law to the data in the 100 MeV-60 GeV 
energy range using a maximum likelihood fit. The dashed blue line presents the exponential cutoff 
power-law model when it is favoured with respect to a simple power-law (TS cuto s > 9, see sec- 
tion (5T)). The statistical errors are shown in black, while the red lines take into account both the 
statistical and systematic errors as discussed in section [BTTl A 95 % C.L. upper limit is computed 
when the statistical significance is lower than 3 a. 
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Fig. 2. — Spectral energy distributions of the off-pulse emission of J2021+4026 (top left), 
J2055+2539 (top right) and J2124— 3358 (bottom), renormalized to the total phase interval. Same 
conventions as for Figure [Q 
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Fig. 3. — Spectral energy distributions of the off-pulse emission of PSR J1023— 5746. The 
LAT spectral points (red) are obtained using the maximum likelihood method described in sec- 
tion [5T27T] in 7 logarithmically-spaced energy bins. A 95 % C.L. upper limit is computed when 
the sta tistical significance is lower than 3 a. The blue points represent the H.E. S.S. spectral 



point s (lAharonian et al.ll2007|) . The Suzaku upper limit is shown with a green arrow (IFujita et al 



20091) . The black line denotes the total synchrotron and Compton emission from the nebula as 



described in section 15.2.21 Thin curves indicate the Compton components from scattering on the 
CMB (long-dashed), IR (medium-dashed), and stellar (dotted) photons. 
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Fig. 4. — Light curves obtained with photons above 100 MeV in a region of 1° around J0034— 0534 
(top left), J0633+1746 (top right), J1813- 1246 (bottom left) and J1836+5925 (bottom right). The 
dashed horizontal line represents the estimated background level, as derived from the model used in 
the spectral fitting. The two dashed vertical lines represent the definition of the off -pulse window, 
as defined in Tabled] Two rotations are shown and 25 bins per rotation. 
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Fig . 5 . — Light curves obtained with photons above 1 00 MeV in a region of 1 around J202 1 +4026 
(top left), J2055+2539 (top right) and J2124-3358 (bottom). Same conventions as for FigureH 
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Fig. 6. — Spin-down flux at Earth as a function of Age for pulsars in the ATNF catalog; Fermi-LAT 
detected pulsars are marked with black circles. Pulsar wind nebulae candidates are marked with red 
stars; LAT pulsars showing a significant off-pulse emission with a plausible magneto spheric origin 
are marked with blue squares. For pulsars with a distance range in Tabled! we use the geometric 
mean of the minimum and maximum values. Note that inaccurate distance estimates can introduce 
artificially low spin-down fluxes, which might account for the handful of pulsar detections below 
10 33 erg s _1 kpc~ 2 . 
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Fig. 7. — Dependence of the PWN luminosity and the pulsar spin down power for LAT detected 
pulsars. Pulsar wind nebulae candidates are marked with red stars; pulsars showing a significant 
off-pulse emission with a plausible magneto spheric origin are marked with blue squares. Only 
pulsars with an estimated distances reported in Table \T\ are plotted. Error bars take into account 
both the statistical uncertainties on the luminosity and the uncertainty on the distance of the pulsar. 
Lines correspond to constant 7-ray efficiency: 100% (dashed), 10% (dotted), 1% (dotted-dashed). 
Pulsars with 2 distance estimates have two markers c onnected w i th gree n dashed error bars. The 
luminosity of the PWN in MSH 15-52 is taken from kbdo etaP (boiOhh . 



